Bouncing transient currents and SQUID-like voltage in nano devices at half filling 



o 
o 

as: 
< 



Michele Cini, 1,2 Enrico Pcrfctto, 3 Chiara Ciccarelli, 1 Gianluca Stefanucci, 1,2 ' 4 and Stcfano Bcllucci 2 

1 Dipartimento di Fisica, Universita di Roma Tor Vergata, 
Via della Ricerca Scientifica 1, 00133 Rome, Italy 
2 Istituto Nazionale di Fisica Nucleare, Laboratori Nazionali di Frascati, Via E. Fermi 4-0, 00044 Frascati, Italy 

Unita CNISM, Universita di Roma Tor Vergata, 
Via della Ricerca Scientifica 1, 00133 Rome, Italy 
^European Theoretical Spectroscopy Facility (ETSF) 

Nanorings asymmetrically connected to wires show different kinds of quantum interference phe- 
nomena under sudden excitations and in steady current conditions. Here we contrast the transient 
current caused by an abrupt bias to the magnetic effects at constant current. A repulsive impurity 
can cause charge build-up in one of the arms and reverse current spikes. Moreover, it can cause 
transitions from laminar current flow to vortices, and also change the chirality of the vortex. The 
magnetic behavior of these devices is also very peculiar. Those nano-circuits which consist of an 
odd number of atoms behave in a fundamentally different manner compared to those which consist 
of an even number of atoms. The circuits having an odd number of sites connected to long enough 
symmetric wires are diamagnetic; they display half-fluxon periodicity induced by many-body sym- 
metry even in the absence of electron-phonon and electron-electron interactions. In principle one 
can operate a new kind of quantum interference device without superconductors. Since there is no 
gap and no critical temperature, one predicts qualitatively the same behavior at and above room 
temperature, although with a reduced current. The circuits with even site numbers, on the other 
hand, are paramagnetic. 

PACS numbers: 72.10.Bg,85.25.Dq,74.50.+r 
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I. INTRODUCTION 



The last decade has witnessed a growing interest in 
the persistent currents in quantum rings threaded by a 
magnetic flux <f>. This problem has many variants. The 
rings of interest may contain impurities, may interact 
with quantum dots or reservoirs. Moreover, both contin- 
uous and discrete formulations have been used to date, 
with similar results. Aligiai modeled the persistent cur- 
rents in a ring with an embedded quantum dot. The- 
oretical approaches to one-dimensional and quasi-one- 
dimcnsional quantum rings with a few electrons are re- 
viewed by S.Viefers and coworkers^; see also the recent 
review by S. Maiti^ where it is pointed out that central 
issues like the diamagnetic or the paramagnetic sign of 
the low-field currents of isolated rings are still unsettled. 
The present paper, instead, is devoted to rings connected 
to circuits. It is clear that the connection to wires sub- 
stantially modifies the problem, but here we point out 
that some of the modifications are not obvious and lead 
to quite interesting and novel consequences. 

Simple quantum rings can be connected to biased wires 
in such a way that the current flows through incquivalent 
paths. In the present, exploratory paper, we consider 
a tight-binding ring with N sites attached to two one- 
dimensional leads and specialize on the half-filled system. 
We show that this innocent-looking situation produces 
peculiar phenomena like bouncing transients i.e. current 
spikes in the reverse direction, and a new sort of simu- 
lated pairing. These systems could find applications in 
spintronic or fast electronic devices exploiting spikes in 
the onset currents, or also in steady current conditions 



when used to measure local magnetic fields. Thus, our 
motivation is twofold. On one hand, we look for the con- 
ditions that can produce novel phenomena in the tran- 
sient current when the system is biased. On the other 
hand we are interested in the bias that develops across 
the system in steady current conditions when it is used 
like a SQUID. After presenting the model in the next 
Section, in Section Mil we present the formalism and in 
Section HV1 apply it to situations where transient currents 
through the circuit are large and appear to bounce to a 
direction contrary to the main stream. We then discuss 
the currents inside the device, with vortices that can be 
clockwise or counterclockwise depending on the parame- 
ters. The nontrivial topology allows interference effects 
of the Aharonov-Bohm type in one-body experiments, 
but in addition here, a many-body symmetry comes into 
play when the ring has an odd number of atoms. The 
interplay of symmetry and topology may lead to a dia- 
magnetic behavior with a half-fluxon periodicity which 
looks like a typical superconducting pattern as shown in 
Section [V] The possible operation of the circuit as a 
magnetometer is also suggested. Our main results are 
summarized in the concluding Section [VTl 



II. MODEL 

The Hamiltonian describing the left (L) and right (R) 
one-dimensional leads is 



H a = t h ^ ( c L,a c m+l,o! + h- c -); 



(1) 
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FIG. 1: Sketch of the site numbering and site energies of the 
triangular ring, shown with part of the leads. The impurity 
is denoted by its energy e c while we shall set e a = e& = 0. 
The hopping integrals in the absence of magnetic field are 

tab — tbc — tea — tdev- 



The electron number current operator between sites 
to and n connected by a bond with hopping integral 
T mn (r m n = *h or t dev e la ^ m ^ inside the device) is 
determined^ by imposing the continuity equation 

Jmn — ~^T~mnc\iCrn ~\~ h.C (7) 

For ring sites, c operators will be replaced by d ones. 
III. FORMALISM 



where th is the hopping integral between nearest neighbor 
sites and c m + i. a annihilates an electron at site to + 1 in 
wire a = L, R. The spin indices are not shown in order 
to simplify the notation. The energy window of both L 
and R continua is (— 2\th\, 2\th\) and the half-filled system 
corresponds to having a chemical potential fj. = 0. In this 
work we consider A^-sided polygonal devices described by 
the Hamiltonian 



JY 



H 



ring 



= ^ (tmndt n d n + h.C.) 



(2) 



<m,n> — 1 



where d n arc Fermion annihilation operators inside the 
device and t mn is the hopping between site to and site 
n. In the rest of the paper we specialize to the case 
tmn = if to and n are not nearest neighbors and t mn = 
tdevG latym,n \ where e IQ ( m ' n ) is a phase factor, otherwise. 
The simplest example is the triangular model [N = 3) 
shown in Figure 1. In this case we also label the sites 
with letters and consider a non-magnetic impurity at the 
vertex site c of strength e c : 



H, 



<m,n>£{a,6,c} 



nnd m d n + h.c.) + e c n c 



(3) 



with n c = d\d c . 

The wires will be attached at a couple of sites in the 
ring that will be specified case by case below by using a 
tunneling Hamiltonian Ht with hopping parameter th- 
Thus, the equilibrium Hamiltonian reads 



Hq — Hl + Hr + firing + Ht 

while for times t > 



(4) 



In the partition-free approach^ the formula for the 
time-dependent averaged current through the m—n bond 
reads 



(J mn (t)) =Tr\pW(t)J mn U(t) 



with 



/(0) 



1 



(8) 



(9) 



the Fermi function computed at the equilibrium Hamil- 
tonian Hq and U(t) the evolution operator. In the actual 
calculations we have adopted a local view and write the 
electron number current 



2e 



{J mil (^)) 

m J2(n\UHt)\r){s\U(t)\m)f(r,s) 



(10) 



where f(r,s) = (r\f^ \s) is the matrix element of /^°^ 
between one-particles states localized at sites r and s. In 
this way we observed that far sites in the wires come into 
play one after another with a clear-cut delay. One can 
simulate infinite leads with wires consisting of I sites and 
obtain quite accurate currents for times up to t = 
the absence of more distant sites does not change the 
results in any appreciable measure. Eventually, when t 
is increased at fixed length /, the information that the 
wires arc finite arrives quite suddenly and a fast drop of 
the current takes placed 

Another useful expression for the bond-current which 
holds for step-function switching of the bias can be ob- 
tained by inserting into Equation (fTDj) complete sets 
■01 , -02 of H eigenstates (in the presence of the bias) with 
energy eigenvalues E(ipi), E(ip2)', one gets 



H{t)=H +H hias (t), 



where 



(5) 



ffbia S (i)= J2 ^ V »(-t)n m ,a, (6) 

a—L.R m— 

with n ma = c\ na c ma . The bias shifts the energy of the 
a sites by a constant amount, with a time dependence 
V a (t). 



2e 



(J mn (t)} = --r Im t mn ^2 e 



(E(i/.i)-E(^2)* 



(n|^)(0 2 |TO) (i>i\ni)(n 2 \i>2}f(ni,n 2 ) ■ (11) 
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This informs us about the frequency spectrum of the cur- 
rent response. The frequencies arise from energy differ- 
ences between the eigenstates of the Hamiltonian of the 
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full circuit with the bias included. The weights at a given 
bond depend in a simple way on the cigcnfunctions at 
the bond and on the equilibrium occupation of single- 
electron states. If there are sharp discrete states outside 
the continuum, (J mn (i)) has an oscillatory component, 
otherwise it tends^ to the current-voltage characteristics 
asymptotically as t — » oo. 



IV. SWITCH-ON CURRENTS 

We study the triangular model of Figure 1 with t a b = 
tbc = tea = tdev and tdev = th for the sake of defmitc- 
ness. Our main idea in the model calculations has been 
to study the e c dependence of the transient in order to 
look for marked out-of-equilibrium behavior, like vortex 
formation, or large charge build-up followed by strong 
current spikes. Also, one is interested in conditions that 
produce a fast change of the response with e c , since pos- 
sibly in such situations magnetic impurities can produce 
strong spin polarization. 

Below we present numerical results for 150 sites in the 
leads (which guarantee an accurate propagation up to 
times 20H/th) and switch on a constant bias Vl — V = 
0.5 and Vr = at t = 0. In the steady state it tends to 
produce a negative local current flow (in the sense that 
the electron number current goes from the right to the 
left wire). 



A. Total current 

In Figure 2 we show the time-dependence of the num- 
ber current in % units on the bond between sites -3 and 
h 

-4 for different values of e c . 



1. Long-time limit 

The times we are considering are short enough to en- 
able us to use the finite-lead scheme specified above, how- 
ever the system already clearly approaches the steady 
state for all e c . The effect of the impurity at site c on 
the asymptotic current is much stronger than one could 
have expected from any classical analogue. The current 
at e c = — 1 (in units of th) is more than an order of mag- 
nitude larger than it is at e c — +1. Moreover, in a range 
around e c = the negative current increases (that is, its 
absolute value decreases) with increasing e c , as one could 
expect if the repulsive site were an obstruction to the 
current flow. However this view is not in line with the 
fact that a further increase of e c increases the conductiv- 
ity of the device. The correct interpretation is that the 
conductivity is ruled by quantum interference between 
the a — b and a — c — b paths, particularly by electrons 
around the Fermi-level. 



2. Short and intermediate times 

Here we are in position to see how this quantum in- 
terference develops in time. It takes a time ~ — for the 
current to go from zero to the final order-of-magnitudc; 
then for some p- large oscillations occur, which appear to 

<-h 

be damped with a characteristic time ~ 10-£. The oscil- 
lations have characteristic frequencies that in accordance 
with Equation (jTTJ) increase by increasing the hopping 
matrix elements or modifying the device in any way that 
enhances the energy level differences. It can be seen that 
the damping of the oscillations is not exponential, and 
characteristic times of the order of 10 ^ are also notice- 
able. 

Figure [2] shows that at e c = the current quickly ap- 
proaches —0.08-^-, but at short times the spread of values 
of the currents is larger than it is at asymptotic times. 
A negative impurity like e c = — 1 produces a spike whose 
magnitude exceeds by ~ 30 per cent the steady state 
value ~ —0.12-^, while a positive e c reduces the conduc- 
tivity of the device. The dependence of the current on 
e c and time is involved, with the e c — 1 case which does 
not belong to the region delimited by e c = and e c = 2 
curves. The e c = 0.8 and e c = 1.0 curves that produce 
small negative currents at long times are most interest- 
ing. They even go positive for some time interval ~ j-. 
Positive currents go backwards. This bouncing current 
is a quantum interference effect, which takes the system 
temporarily but dramatically out of equilibrium and in 
counter-trend to the steady state. The main spike of re- 
versed current is much larger than the long-time direct 
current response. Next, to understand what produces 
the bouncing current, we look at the transport inside the 
triangular device. 
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FIG. 2: Time-dependence of the current on the bond between 
the sites -3 and -4. The numbers shown close to each curve 
indicate the values of e c . From bottom e c = —1,0,2,08, 1 in 

units of th- Here and below t is in -J*- units and the current 

h 

is in units of -jf . 
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FIG. 3: Time-dependence of {J a b}- The numbers shown close 
to each curve indicate the values of e c . From bottom, e c = 
2, -1,0, 1, 0.8 in units of t h . 



porary back- flow. After the burst, J a b remains positive 
(e c = 0.8 and e c = 1.0 curves in Figure 3) while J c (, re- 
mains negative value (e c = 0.8 and e c = 1.0 curves in 
Figure 4). In other terms, we observe the formation of 
an anti-clock-wise current vortex. 

Remarkably and unexpectedly, by increasing e c one 
reaches a critical value beyond which the vortex becomes 
clock-wise, as one can see from the e c = 2 curves in Fig- 
ures 3 and 4. This conclusion is unavoidable since the 
current changes sign in both arms of the circuit. The to- 
tal current (Figure 2) remains negative, as one expects. 
We were unable to find any simple qualitative explana- 
tion for the inversion of the vortex. In general the criti- 
cal value of e c depends on the bias V and for V = 0.5 is 
about 1. Moreover, we emphasize that the current across 
a strongly repulsive site with e c = 2 is still comparable 
in magnitude with the one on the a-b bond. 




B. Laminar flow and Vortices 

In Figures [3] and |4] we report the time-dependent {J a b) 
and (Jbc) at times < t < ^ for e a = = and various 
values of e c . When the site c is attractive for electrons 
(that is, e c < 0), despite some transient oscillations, the 
current remains negative, both on the a-b and c-b bonds; 
the current flows from b towards a and c, so the flow 
is laminar. The magnitude of the current on the a-b 
and c-b bonds is comparable. Although site c may have 
a high electron population, the local current does not 
concentrate on either bond. 

For e c = 0.8 and e c = 1.0 the current J a b on the a-b 
bond (Figure 3), after a negative transient spike, pro- 
duces a positive one, which is the main contribution to 
the bouncing current noted above. This behavior can be 
understood in terms of a strong charge build-up on the 
a-c-b arm of the ring during the first burst following the 
switching of the bias, which eventually triggers the tem- 



FIG. 5: The -| periodicity arising from the symmetry-induced 
false pairing. Top left: the triangular device connected to 
A =14-site wires. The other panels show the ground state 
energy Eo(a) versus a. Top right: tdev = th- Bottom left: 
tdev = 3t h . Bottom right: t dev = 5t h . 

Nanoring devices asymmetrically connected to wires 
of A sites each are even more peculiar for their magnetic 
properties. Here we restrict to the case e c = 0. The 
magnetic flux <p through the ring is inserted by the Peierls 
prescription™ 

where a = and (f>o = is the magnetic flux quan- 
tum. In the case of vanishing bias, in stationary condi- 
tions one finds a diamagnctic current which is confined 
to the triangular ring. In Figure 5 we illustrate our re- 
sults. We sketch a triangular device connected to 14-site 
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leads and the ground state energy E${a) versus a = 4- 
for tdev = tab = the = t ca = 1,3 and 5. The expected pe- 
riodicity in a with period 1 due to gauge invariance is of 
course observed. Strikingly however, we notice that the 
period is actually ^, and the system is diamagnctic, i.e. 
the energy increases when the un-magnetized system is 
put in a magnetic field. In other terms, Er,(a) looks like 
the ground state energy of a superconducting ring, al- 
though this is a non-interacting model and the spectrum 
is gapless. 

One could wander how the positive diamagnctic re- 
sponse at small fields arises, since the initial dependence 
of the total energy on the flux is quadratic and the 
second-order correction is always negative. However this 
is an apparent paradox. The energy change is due to a 
perturbation td e v{e 2ma — 1) ~ t<i ev (2nia — 2ir 2 a 2 ) which 
is the effect of changing the phase of one bond. First or- 
der perturbation theory corresponds to compute the cur- 
rent operator over the ground state and it is zero. The 
quadratic contribution is given by second-order pertur- 
bation theory in 2matdev plus a term coming from the 
first-order correction in — 2i: 2 a 2 tdev The former term 
is always negative while the second term can be cither 
positive or negative. 



A. Nontrivial role of the wires 

The results in Figure [5] are striking because the iso- 
lated nano-ring with 3 sites at half filling does not sim- 
ulate any superconducting behavior; instead, it yields a 
paramagnetic pattern, symmetric around a = ^ (a = 1 
is equivalent to a = 0). One can easily work out the 
lowest energy eigenvalue Eo with 3 electrons. Eq = — 3 
for a = sinks to Eo = — 2y3 ~ —3.46 at a = and 
raises again to Eq = —3 at a = \. Thus, there is a half- 
fluxon periodicity at half filling, but a = and a = h 
are maxima and correspond to degenerate 3-body states. 




FIG. 6: Dependence of the barrier hight B on the total num- 
ber N of sites, which include the triangular device and the 
L and R leads. The parameters are the same as in Fig. 6. 
The line is a guide for the eye. Note that for short leads B is 
negative. The inset shows the almost linear dependence of B 

On tdev — t a b — t a c — the- 

The supcrconductor-likc response requires the presence 
of wires, despite the fact that the diamagnctic currents 



induced by the field are strictly confined to the triangu- 
lar device. The currents do not visit the wires, but the 
electron wave functions do. In order to produce the dou- 
ble minimum, the wires must be rather long. The barrier 
height B = Eq(j) — Eq(0) depends on the total number 
of atoms Af = 3 + 2A (see Figure 6) and below a mini- 
mum length B < 0. In Fig. [5] we plot the dependence of 
the barrier B on the number of sites Af = 3 + 2A of the 
leads. We observe that B saturates with increasing Af. 
This finding is noteworthy and unusual. The inset of Fig- 
ure 6 shows how B depends on \td ev \ and suggests that 
except for an initial quadratic region the dependence is 
basically linear. With hopping integrals in the eV range 
B easily exceeds room temperature. 



B. Bipartite and not bipartite wired devices 

Next we investigate why the effect takes place in this 
geometry and at half filling. A crucial observation is that 
the system depicted in Figure 5 is not a bipartite graph. 
By contrast, Figure 7 shows the flux dependence of the 
ground state of a bipartite graph, namely, a square device 
connected to wires. The response is paramagnetic since 
the ground state at cf) = is degenerate, the field lifts the 
degeneracy, a Zeeman effect occurs, and the ground state 
energy is lowered. In addition, the trivial periodicity is 
observed. In Figure 8 we add one site to the device, 
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FIG. 7: Flux- dependent energy Eo(a) for a square device with 
tdev = th- Notice the scale, the paramagnetic behavior, and 
the level crossing at no flux. The central part of the system, 
with the device and the beginning of the wires, is also shown. 

and the diamagnetic double minimum pattern is found 
for the resulting pentagonal device, which does not pro- 
duce a bipartite graph. Based on this observation and on 
a symmetry analysis, we can show that the half-fluxon 
periodicity {E {\) = E Q {Q)) holds. This is a theorem 
which holds for any ring with an odd number of atoms 
connected to leads of any length provided that the site 
energies vanish. 



C. Symmetry analysis 

Let C denote the charge conjugation operation (or 
electron- hole canonical transformation) c a — > b' a , where 
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FIG. 8: Flux-dependent energy Eo(a) for a pentagonal device 
with tdev = th- Notice the diamagnetic behavior, and the 
absence of level crossing at no flux. The central part of the 
system, with the device and the beginning of the wires, is also 
shown. 

c a annihilates electrons and creates a hole with quan- 
tum numbers a. C is equivalent to r„ m — > — r mn 
throughout, or, since the site energies vanish, to H — > 
-H. 

We recall two elementary results of graph theory: 1) 
an isolated ring is a bipartite graph if and only if it has 
an even number of atoms 2) adding wires of any length 
does not change the result. 

In bipartite lattices, C is equivalent to a sign change 
of alternating orbitals, which is a gauge transformation. 
Hence, H and —H have the same one-body spectrum, 
that is, the spectrum is top-down symmetric. 

For the wired triangular ring and any other non- 
bipartite graph, the one-body spectrum is not top- 
down symmetric (except^ at a = \,ot = | ) so it 
does not appear the same after the transformation, and 
C : H — ► — H is equivalent to changing the sign of one 
bond in the ring. But this is just the effect of the oper- 
ation F : a — > a + | which inserts half a fluxon in the 
ring. Therefore the combined operation C o F is an exact 
symmetry of the many-electron state which holds at half 
filling. It is clear that F : H — > —H i.e. F turns the 
spectrum upside-down. 

However, as noted above, the spectrum is not top- 
down symmetric and when at </> = ^ the occupied and 
empty spin-orbitals have exchanged places the system 
does not quite look like it was at a = 0. For instance, 
at the top of the spectrum of the device of Figure 5 at 
a = there is an empty split-off state. At a — 0.5 this 
becomes a deep state below the band. Nevertheless, we 
wish to prove that the many-body ground state energy 
Eo(a) at a — | is exactly the same as at a = 0. 

The reason lies in another symmetry of the many-body 
state at half filling. In terms of one-body spin-orbital 
levels, Tr-ff(a) = implies 

E (a) = Ma)e l (a) = - ^(1 - n;(a))e;(a). (13) 

i i 

Since under F the negative of the energy of the unoccu- 
pied levels coincide with the energy of the occupied ones, 
we have Cj(^) = — Cj(0) and (1 — ra;(|)) = nj(0) from 
which it follows that 

£0(5) = -£(1- ^(^ = £0(0). (14) 

i 




FIG. 9: Operation of the triangular device with e a = et = 
e c — and tdev ~ th used a a SQUID thread by a flux cj>. 
A fixed current J = 0.016%- flows through the device. The 
right wire site energies are raised by Vr = V(4>), the left wire 
site energies are lowered by Vl = —V((j)), and V is adjusted 
in order to keep J fixed. The V(4>) versus <j> plot is perfectly 
periodic with a half-fluxon period and simulates a supercon- 
ducting SQUID, although no superconductors are needed. 
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FIG. 10: Here we present the same thought experiment as in 
the previous Figure, except that the device is a square with 
tdev ~ th threaded by a flux <f> with zero on-site energies. The 
SQUID-like behavior and the half-fluxon periodicity are lost, 
as expected in bipartite graphs (see text). In this particular 
example the lower-left bond was raised to 1.2th in order to 
study the effects of a slight deformation of the square, but 
no relevant departure from the perfect square behavior was 
found. 



D. Superconductor- free Quantum Interference 
Device 

The superconducting quantum interference device 
(SQUID) consists of two superconductors separated by 
thin insulating layers to form two parallel Josephson 
junctions, and can be used as a magnetometer to de- 
tect tiny magnetic fields. If a constant biasing current 
is maintained in the SQUID device, the measured volt- 
age oscillates with the change in the magnetic flux, and 
by counting the oscillations one can evaluate the flux 
change which has occurred. In principle one can produce 
Josephson-like oscillations without the need of supercon- 
ductivity; if the device can be realized this can be a phys- 
ical principle of interest for applications, maybe even at 
room temperature. Impurities and any disturbance af- 
fecting the quantum coherence can lower the operating 
temperature, but since in principle B can be as large as 
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leV we believe there is enough motivation for an exper- 
imental activity on this idea. According to the above 
arguments, the triangular pentagonal and other odd- 
numbcrcd polygonal devices should simulate a SQUID, 
while bipartite graphs should present a normal behav- 
ior. We have performed the thought experiment with 
the triangular device connected to infinite wires; the re- 
sults arc reported in Figure [5J We keep a fixed current J 
flowing through the device by adjusting V(<p), while the 
flux (j> threading the device is varied. It can be seen that 
the plot of V(4>) versus <f> is periodic with a half-fluxon 
period. The system simulates a SQUID, although no su- 
perconductors are needed. A counter-example is given in 
Figure [101 where the triangular device is replaced by a 
square one and the effect disappears. 

VI. CONCLUSIONS 

We have shown that simple asymmetric closed circuits 
have rather subtle properties when unsymmetrically con- 
nected to the biased circuit, that could also be useful for 
designing new kinds of devices. We have found that an 
impurity site of energy e c in the longer arm of a triangu- 
lar device can cause a transient bouncing current, which 
goes in the opposite direction than the long-time current 
and is much more intense but lasts for a short time. Look- 
ing at the current distribution inside the device, e c < 
favors a laminar current flow; instead, e c > produces 
an anticlockwise current vortex; when a critical value is 
exceeded, however, the vortex chirality reverses. Fur- 
ther, we have shown that a class of closed circuits at half 
filling have a degenerate ground state which is paramag- 
netic, i.e. gains energy in a magnetic field by a Zeeman 



splitting. The magnetic behavior is completely changed 
when the circuits are connected to leads and form a non- 
bipartite lattice. Long leads constitute an essential re- 
quirement and although the diamagnetic currents are 
confined to the ring, the leads modify the magnetic prop- 
erties substantially. With long enough wires one obtains 
a diamagnetic behavior with half fluxon periodicity and 
a robust barrier separating the energy minima at and 
4r . This pattern mimics a superconducting ring although 
there are no gap, no interactions, no critical temperature. 
We have traced back the origin of this fake superconduct- 
ing behavior to a combination of charge-conjugation and 
flux which provides a symmetry of the many-electron de- 
tcrminantal state. Finally we have shown that in princi- 
ple one can extend this simulation to the point of build- 
ing a functioning interference device capable of measur- 
ing local fields and analogous to a SQUID but working 
even above room temperature. It cannot be excluded 
that even-sided circuits can be useful for the same pur- 
pose, and indeed Figure 10 suggests that one could do so, 
exploiting the trivial periodicity. The signal in Figure 9, 
however, is much more monochromatic, and this suggests 
that the odd-sided version should make it much easier to 
read-off the magnetic field intensity from the amplitude 
of the voltage oscillation. 
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